Ferromagnetism and superstructure in Cai.^La^Bg 
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We critically investigate the model of a doped excitonic insulator, which recently has been 
invoked to explain some experimental properties of the ferromagnetic state in Cai-^La^Be. We 
demonstrate that the ground state of this model is intrinsically unstable towards the appearance of 
a superstructure. In addition, the model would lead to a phase separation regime and the domain 
structure which may be prevented by the Coulomb forces only. Recent experiments indicate that a 
superstructure may indeed show up in this material. 

PACS numbers: 71.10.Ca, 71.35.-y, 75.10.Lp 
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The discovery of weak ferromagnetism in lightly-doped 
hexaborides Cai_ x La x B6 [lj] has renewed theoretical in- 
terest in the so-called "excitonic" transition vigorously 
discussed in the 60-s and 70-s g-g](see for a review). 
Band structure calculations ||] indicate that hexaborides, 
DBg, may, in fact, be semimetals owing to an acciden- 
tal small band overlap at the three X-points of the Bril- 
louin Zone, two bands at each X-point having symme- 
try X 3 , X' 3 . It was suggested |J that the use of the 
model 0] together with the energy mechanism for ferro- 
magnetism first suggested in Ref . [10 11 provide at least 
a qualitative explanation for the surprising findings EJ. 
Bearing much of the resemblance with the mathematics 
of the BCS weak-coupling theory of superconductivity, 
the Keldysh-Kopaev model is a convenient tool which 
makes possible to treat an excitonic transition in the con- 
trollable manner. 

Assuming that the basic physics of hexaborides is prop- 



erly accounted for by this oversimplified weak coupling 
scheme, in what follows we study in some more details 
the zero temperature phase diagram as a function of the 
doping level, and its' stability to anisotropy features. The 
main result is that the system almost inevitably devel- 
ops a superstructure on the background of the initially 
cubic lattice. Our explanation for the magnetic moment 
per doped lanthanum ion to vary and even become very 
small differs from the one suggested in Ref. |g|. 

Zhitomirsky et al. [|J have used the model of Ref. || 
of an excitonic phase transition in a semimetal to explain 
weak ferromagnetism observed in hexaborides. Ferro- 
magnetism is known to appear in this model, and was 
investigated in detail by Volkov et al. [J10],[ll]]. We first 
briefly recall the mechanism of magnetic moment forma- 
tion in this model. It is well known 0] that if the screened 
Coulomb interaction between an electron and a hole in 
two bands, 
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is dominant, and the electron and hole Fermi surfaces 
do coinside ("nesting"), the excitonic transition has 
a degeneracy for the onset of CDW ("singlet") and 
SDW( "triplet") excitonic condensate. This degeneracy 
is lifted only by additional weak short-range Coulomb 
terms, which favor triplet excitonic condensate pip), or 
electron-phonon interactions, which favor singlet (CDW) 
state. This splitting, however, being usually considered 
to be weak at 5g ~ g 2 , where g — V(0)N n <C 1 is the 
screened Coulomb coupling constant, N — rnkp/i^ 2 ) 
is the density of states per spin for a single band, the tem- 
peratures of the triplet or singlet excitonic transitions 
are close (on the exponential scale, T c ~ exp[— ff -1 ]), 
T s o ~ T t o- Then, as it was first shown in Ref. [pOlJTTj and 
applied to the physics of hexaborides in Ref. [BL ferro- 
magnetism can appear as a result of doping due to the 



development of a triplet excitonic instability in the pres- 
ence of a singlet order, or vice versa. Ferromagnetism is 
the direct result of the fact that in the presense of these 
two condensates both crystalline and time-reversal sym- 
metries get broken. Summation of the leading logarith- 
mically divergent terms in the presence of one condensate 
explicitly confirms the divergent Curie- Weiss behavior 
of the spin susceptibility with some Curie temperature 
found in Ref. 0. 

For the analysis of ferromagnetism in the doped state 
at low temperatures one may neglect the small differ- 
ences in the energy spectrum of the SDW or CDW ground 
state B . The energetic analysis is straightforward in case 
when triplet and singlet coupling constants are equal M , 
since then equations for different spin polarizations are 
decoupled. Indeed, for a simple model with two isotropic 
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and e e (k) = |- - /j, 



1m f* ~ 2 ' 

., i f(+ ^, the zero temperature excitonic 
gap in the unpolarized state is given by 



e. 



A 2 = A (A -2n), 



(2) 



where A = 2e c exp(— g^ 1 ) is the excitonic gap at zero 
doping, e c is a cutoff energy around the Fermi surface, 
n = N/ANq is the concentration of doped carriers in en- 
ergy units (the level of the chemical potential in the bare 
metallic phase is fj, = n). The ground state energy per 
spin direction relative to the stoichiometric normal state 
is given by: 



1 



E(n) = N n' ~ -N (A - 2n) z . 



(3) 



Thus, the excitonic pairing disappears at n = 0.5Ao. A 
quick analysis || shows that the energy of a polarized 
excitonic state, E[n(l + p)] + E[n(l — p)] has a minimum 
value for p = 1, i.e. for the complete polarization of 
added carriers. Therefore, below T c , there are two order 
parameters, A s and A t , and, correspondingly, two differ- 
ent gaps in the spectrum for different spin polarizations: 



v / A (A -4n) < n < A /4 
A T = A < n < A /2. 



(4) 



Hence at T=0 electrons and holes are paired for only one 
spin direction when Ao/4 < n < Ao/2. (The system 
undergoes a 1-st order phase transition into the unpolar- 
ized normal metal at n cr = Ao/2.) Thus, this mechanism 
gives a large effective moment equal to \^,b per doped La 
atom, and some efforts have been applied in Ref. |1 to 
argue that this moment may be forced to become small if 
the interaction between the orientation of the magnetic 
moment and the direction d of the triplet order param- 
eter (SDW), A t (p) = (d(p)o - ), was taken into account. 
The mechanism of Ref. || seem to us too artificial, first 
of all, because the energy difference between CDW and 
SDW states was ignored. Meanwhile, it is easy to see 
that an anisotropy of electron and hole pockets would 
reduce the net magnetization. Indeed, the two oppo- 
site spin polarizations preferred by the system are gov- 
erned by the position of the chemical potential. First, an 
anisotropic solution for the order parameter, A(p), would 
result in variation of the gap itself along the Fermi sur- 
face. Secondly, even a small anisotropy ( "antincsting" ) of 
the electron and hole spectra hinders the excitonic gap, 
even leading to the gapless pockets along the initial Fermi 
surface [0,13[. Therefore, while starting from the spin- 
up and spin-down spectra for the isotropic gap shown 
in Fig. la (1 /ib per La), one would end up with doped 
electrons to spill over between the two energy branches 
thus obviously reducing the total magnetization, as it is 
shown in Fig. lb. 
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FIG. 1. a) Energy bands and spin directions in the isotropic 
excitonic insulator; b) the value of the total magnetic moment 
depends on the shape of the energy spectrum in momentum 
space along the Fermi surface when the gap anisotropy and/or 
"antinesting" are taken into account. 

We will not pursue the detailed calculations for these 
mechanisms, since the weak coupling model |3| of an exci- 
tonic transition suffers from several obvious defneiencies, 
each of which results in an instability towards a forma- 
tion of an inhomogeneous structure. First of all, for the 
homogeneous excitonic state to exist, the electron and 
hole Fermi surfaces should be sufficiently close to nest- 
ing. Indeed, a quick calculation shows that in a model 
with mass anisotropy, 
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a homogeneous excitonic phase disappears at 
5m 2eAg 
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where Ao is the excitonic gap for the isotropic situation 
(5m = 0); e = 2.71828. Since Aq is exponentially small 
in terms of e c ~ E gi Aq — 2e c exp(—l/g), the weak cou- 
pling model allows only a modest mass anisotropy. 

The range of anisotropy for exciton formation may 
be extended beyond Eq.(pf), however, since at larger 
anisotropies it leads to an inhomogeneous state with long- 
wavelength oscillations of spin and (or) electron densities 
|14| . Similar drawback of the Keldysh-Kopaev model be- 
comes obvious when one considers doping even in a model 
with the perfect nesting. Mathematically the problem in 
the last case becomes equivalent to the problem of the co- 
existence of the superconductivity and ferromagnetism 
|15[ . The exchange field, /, in that model is equiva- 
lent to the chemical potential fj, in the excitonic state. 
Thus, the gap equation has two solutions: 1) A = Ao; 
2) A 2 = 2/A - A 2 ,, where A is the gap at / = 0. 
The energies of the ground state are easily found: 1) 
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The branch (2) is not stable for the homogeneous su- 
perconductor iQ. The physical difference between the 
problems of the doped excitonic insulator and the ferro- 
magnetic superconductor is that in case of an excitonic 
insulator the number of added dopants is fixed, not the 
chemical potential. Then the "stable" solution (1) cor- 
responds to zero doping, the case when the chemical po- 
tential fi lies below the gap edge A . The dependence 



J-c = exp 



of a homogeneous solution for the gap on doping then is 
given by (2), with jj, being expressed in terms of added 
dopants, n = y/i 2 — A 2 (where n is in energy units), re- 
produces Eq.(g) and the effects considered by Volkov et 
al. [|l(j,[ll| and Zhitomirsky et al. Q. 

The instability of the homogeneous solution at larger n 
is seen from the form of T c (n), explicitly calculated (see 
Fig. 2) in Ref. f§: 
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At high enough densities, n > A /2, T c (n) displays a 
reentrant behavior from the excitonic insulator into the 
metallic state at low temperatures. The behavior indi- 
cates some sort of an instability or a 1-st order phase 
transition. However, similar to the case of non-perfect 
nesting above, one can again check for an instability to- 
wards a transition into an excitonic state with an in- 
commensurate wave vector q, searching for a maximum 

T c (|q|). 
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FIG. 2. Transition temperature into the excitonic state as a 
function of the concentration of dopants, n = iV/(47Vo). The 
dashed line shows the phase transition into inhomogeneous 
excitonic insulator with the wave vector q. 

The result of our numerical calculation is shown in 
Fig. 2. At low dopings T c corresponds to the homogeneous 
state with q = 0. The dashed line in Fig. 2 shows the 
temperature for the onset of an inhomogeneous phase, 
T c (|q|), where the value of |q| itself depends on n. (At 
low T c |q| = 2An/vp.) The excitonic regime at T = 
for such inhomogeneous state appears by a first order 
transition at n* — 0.71A (see below) and extends to 
n cl ~ 0.755A . (Inhomogeneous state was first discussed 
by Rice in connection with itinerant antiferromagnetism 
in chromium fl7|| .) 

Assuming again that CDW- and SDW- states have 
equal energies (Sg = 0), it is easy to show that T c in 
the presence of magnetic field, B, initially increases with 



B, as follows from T c (n, B) obtainable by a mere substi- 
tution: 



T c (n, B) = T c (n - /z B B), 



(8) 



from Eq.(0). The |q|-value in the presence of the mag- 
netic field B also follows from the same substitution. 
This result may also be considered as another manifes- 
tation of the system's tendency towards a ferromagnetic 
state. According to Eq.(g), at a higher doping the mag- 
netic field may cause reentrance into the excitonic insu- 
lator phase. 

The vicinity of the upper concentration, n c \ , where the 
non-homogeneous solution first appears if the concentra- 
tion is decreased from the side of the normal metal, may 
be studied in the same manner as for the corresponding 
superconductivity problem, producing the well-known 
Larkin-Ovchinnikov-Fulde-Ferrel state p8J19|. The so- 
lution in our case (we need to find /i(n)) again secures 
the stripe phase as the most energetically favorable one. 
The density of carriers then oscillates according to: 



n(x) = n - n[A 2 /n 2 1 ](l + 1.3cos(2gx)), 



(9) 



where qv = 2.4(n c i + 5n/l5.7), A(x) — 2Ai,cos(qx), and 
|A 6 | 2 = 0.936n c i(n cl -n). (10) 

Here v is the Fermi velocity. 

An inhomogeneous distribution Eq.(J9|) of a weak 
enough charge on the scale of a coherence length £0 > 
^tj?, where !;tf is the Thomas- Fermi screening radius 
should not present a major trouble, since the lattice can 
adjust itself to produce a periodic modulation. Yet it 
remains unclear how at low T the system evolves when 
one proceeds from the metallic end by further decrease 
of the concentration of dopants. 

In any event, if the Coulomb forces are completely 
neglected, with further concentration decrease the ho- 
mogeneous phase is restored through the phase sepa- 
ration regime. This starts to take place at the point 
n* = Ao/\/2 Em. To find the two phases and how they 
coexist, in our case turns out to be completely equiva- 
lent to minimizing the energy of the intermediate state 



of the 7 s *-type superconductor in a fixed magnetic field. 
The same energy considerations show that the system 
separates into the mixture of two phases: the one with 
n e = Uh and the other with n e — Uh = n* = Ao/v2. The 
relative volume fraction of each phase is given by: 



V x /V = 1- n/n* V 2 /V = n/n* 



(11) 



The domain sizes would be determined by the surface 
energy. Without the bulk Coulomb forces such phase 
separation is energetically even more favorable than fer- 
romagnetism. On the other hand, the Coulomb energy 
would work against the spatial charge separation, which, 
hence, may stabilize a homogeneous regime in a range of 
small enough n. In view of these uncertainties we made 
no attempt to address the other issue, namely, that for 
a complete description of the true ground state of DBg 
in the frameworks of the model |j|| , one has to account 
for the three X-points in the cubic system related to each 
other by the symmetry transformations. This has been 
done recently in Ref. fl20|j for the problem of multiband 
superconductors, and the results of Ref. (20| are imme- 
diately mapped onto the problem of excitonic insulator 
with the spectrum of Ref. M . 

In conclusion, we have shown that the excitonic model 
p|, even though may capture a qualitatively correct 
physics, necessarily leads to appearance of a superstruc- 
ture as far as the doping dependence is concerned. We 
are aware of the only experimental observation yet of 
such a superstructure Q which, at least, does not con- 
tradict our expectations above. We emphasize again that 
our explanation for possible small magnetic moment dif- 
fers from the one given in Ref. |1. The magnetic mo- 
ment can be as big as ~ \[Xb per doped La in the model 
JlO| |. Effects of anisotropy and changes in the energy 
spectrum through doping would change this sig nificantly, 
as mentioned above. Recent experiments []22] on BoBq 
doped with La have produced moments about 0.4/is at 
x = 0.05. 

After this work was performed, the authors have dis- 
covered the preprint Ref. [£3| by Balents and Varma. We 
shall not discuss their results regarding the three different 
X-points. We note that our results on the doping depen- 
dence agree for the most part. We do not agree, however, 
with the claim that the results of |H|[ll],|| contain signif- 
icant errors, since the macroscopic phase separation is to 
be made difficult by the Coulomb terms. 
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